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Introduction 



Recent theoretical studies show that the properties of nuclear matter can be described 
nicely in terms of the Relativistic Mean Field Theory (RMF) |[l|, |2|, Guichon proposed 
an interesting model on the change of the nucleon properties in nuclear matter (QMF or 
quark-meson coupling model (QMC)) 0. The model construction mimics the relativistic 
mean field theory, where the scalar a and the vector meson 00 fields couple not with nucleons 
but directly with quarks. The quark mass has to change from its bare mass due to the 
coupling to the a meson. In this work we shall investigate the quark matter within the 
Relativistic Mean Field theory motivated by the Nambu-Jona-Lasinio (NJL) model |[]||[]]|[]- 
The three-flavor NJL model has been discussed by many authors, e.g. |||]. To include the 
strange quark the QMF model is extended by incorporating an additional pair of the strange 
meson fields which couple only to the s quark It is very interesting to construct the 
Guichon model, where the nucleon is described in terms of constituent quarks, which couple 
with mesons and gluons. This model (we refer it as the quark mean field (QMF) model was 
named in ||10|)- The meson fields act on quarks inside a nucleon or the neutron or quark 
star and change the bulk nuclear properties. The QMF model predicts an increasing of the 
nuclear size and a reduction of the nucleon mass in the nuclear environment of the star. 

The properties of strange quark matter (SQM) and the conjecture that SQM could be the 
absolute ground state of strongly interacting matter ||Tl|, |12| have attracted much attention 
in nuclear physics and astrophysics. In 1984 Farhi and Jaffe ||T^ found that SQM is stable 
compared with an ^^Fe-nucleus. Still the bag model calculations play a key role in this field 

0- 

In the framework of the NJL model Babulla and Oertel § have showed that due to a 
large constituent quark the strange star is not absolutely stable. The quark-meson coupling 
will reduce the effective quark mass what may make the existence of the strange star more 
possible. 

The aim of this work is to examine the infiuence of the quark-meson coupling model on 
the quark star or neutron star having the quarks core. 
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The Quark Mean Field Theory 



The fields of the model QMF for a, uo and p-mesons are denoted as cu^, p^. To reproduce 
the observed strongly attractive A A interaction two additional meson fields are added, the 
scalar /o and meson fields denoted as the fields V9*and 0^. The Lagrange density function 
for this model has the following form 

^j:}=iL~frd^.Lf - j:}=i9f(LfHeRf + h.c.) + 
iqYD^q - q{mQ - g^iph - gl^^J^)q+ 

The field stress tensor of the vector mesons have the following structure: 



(1) 



(2) 



(3) 



The covariant derivative acts on baryons as follows: 



(4) 



The potential in RMF ||T|, [T^, is given by 

U{^) = Bo + ^mlif^ + ^g^v^ + ^g3fi 
The fermion fields are composed of quarks and electrons, muons and neutrinos 
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Parameter 


TM1|1.5, 16| 


M 


938 MeV 




783 MeV 


Mp 


770 MeV 




511.198 AfeF 




975 MeV 




1020 Mey 




310.6 MeF 


TTin = m* = m*, 7^ m* 


367.7 MeF,549.5Mey 0] 




7.2325 fm-^ 


54 


0.6183 


5s = 5o-Af = 35^ 


10.0289 


gw = gu)N = ^gZ 


12.6139 


gp = gpN = 5p 


4.6322 


5s* = 5o-*A = '^/'igs 


5s* — 35's* 


50 = 5M = V^/^gw = ^gl 


C3 = 71.3075 



TABLE I: Parameter sets for the Lagrangian function (|l|). 
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nis, Mo;, Mp are masses assigned to the mesons fields Here q denotes a quark field with three 
flavors, M, d and s, and three colors. The Lagrangian function includes also the nonlinear 
term ^Cs{uj^lj^)^ . We start with equal dynamical quark masses niQ = m* = = m* = 
310.6 MeV (the 5'f/(3)-symmetric quark matter with starting mo,u = '^o.d = "^o,s = 0)Q 
and a bag pressure Bq = 57.3 MeV/iin^ originated from the Nambu-Johna-Lasinio model 

Lnjl = iqi^D^q - qm^q + G ^{(gA^g)^ + (zgA^^g)^}. (9) 

a=0 
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We restrict ourselves to the isospin SU(2) symmetric case, ttiq = ttiq, so 



mo 



\ 



V 



mo,s j 

In the quark massless limit the system has a f/(3)L x f/(3)i? chiral symmetry. The model is 
not renormalizeable and we have to specify a regularization scheme for divergent integrals. 
For simplicity we use a sharp cut-off A in 3-momentum space. Besides A we have to fix 
the coupling constants G and the current quark masses mo,^ and mo,s. In most of our 
calculations we will adopt the parameters of g A = 602.3 MeV, GA^ = 1.835. With 
A, G and K as specified above, chiral symmetry is spontaneously broken. The possible 
existence of deconfined quark matter in the interior of neutron stars using the XJL model 
was examined in paper |^. In the work the non-symmetrical quark masses vtlq = {m* = 
= 367.7 MeV, m* = 549.5 MeV} are obtained (with starting mo,u = mo,d = 5.5 MeV, 
mo,s= 140.7 MeV). 

The Euler equation for $yi = {(/?, ip^,, lj^, p^, (j)^, q} fields are 



Dip = mlif + g2(p'^ + gsif^ 



gUhq 



(10) 

(11) 



where 



- d.r'' = Mir - giJun) 



J: 



5C/(2) 



qYhq 



is the SU(2) up and down quark current while 



(12) 
(13) 

(14) 



Ju{i) = qYhq 



(15) 



is the U(l) strange quark current. The nonabelian gauge field p*^'" obeys 



- D.R^'''^ = M' p-'- - gU- 



(16) 



where 



Js = ^qYr'q 



(17) 



is the isospin SU(2) current. In the system we have conservation of "nucleon" charge 



73 _ tO 

'SU{2) 



the isospin charge 

and the "hyperon" charge 



Qh = I d^xJ[ 



3^ 7O 
Uil) 



The last is the Dirac equation 

tYD.q - (mg - gtvh - gl^.h)q = 0. (18) 
The physical system is totally defined by the thermodynamic potential ||l] 



Vt = -fcT/nTr(e~'^(^~''^^'^"^^^3^''«^«)) (19) 
where H is the Hamiltonian of the physical system 

H = J2 [ d'x{do<^A4 - ^} (20) 

and vr"^ = g^^^^-, is a momentum connected to the field The fields $a = y^^,, u;^, 
P^n 0/i7 9} denote all fields in the system. In this paper we shall use the effective potential 



approach build using the Bogolubov inequality |15, 18 



Vt<Q.i = VtQ{mB,mF)+ < H - Hq>q (21) 

VLq is the thermodynamic potential of the trial system as effectively free quasiparticle system 
described by the Lagrange function 

-IKR^"" - W.^F'' - (22) 

Similar to the general case 
and 
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and for gluons 



We decompose the $a field into two components, the eff'ectively free quasiparticle field $^ 
and the classical boson condensate 

$A = $A + U (23) 
In the case of the RMF model we have 

<^ = ^ + a (24) 

V^* = ^* + cr* (25) 

= + ^A" = '^M.ow^ (26) 



Pl=r^rl, r; = 6'^%,or (27) 



4>fi = 4>i, + w^M*, w^M* = ^fifiW* (28) 

The = {o", w^, w^*} field will be treated as the variational parameters in the effective 
potential. Also the boson and fermion mass ttib, mp will be treated as as the variational 
parameters. The covariant derivative for the trial system is 

D^ = d^ + UgyX + ^^>M + (29) 

This introduce the homogenous fermion interaction with boson condensate w^, r^, w^*. The 
fermion quasiparticle will obey the Dirac equation 

{i^^D^-mF)^pf = Q (30) 

The constant condensate w, r simply shift the chemical potential from fii = /i° (when 
w = r = 0) to 



1 

2- 
1 

5 

= - glw^ (33) 



= ^° + 2^p'' ~ 9> (31) 
= /^d - - glw (32) 



where fiu = + |/^3 and fid = fJ' ~ |/U3- 

Quarks and electrons are in /5-equilibrium which can be described as a relation among their 
chemical potentials 

where fid, f^s and fie stand for quarks and electron chemical potentials respectively. If 
the electron Fermi energy is high enough (greater then the muon mass) in the neutron star 
matter muons start to appear as a result of the following reaction 

c? — > n + e~ + I7e 
s — > n + /i" + 17^ 

In a pure quark state the star should to be charge neutral. This gives us an additional 
constraint on the chemical potentials 
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-^riu - -Ud - -Us - rie = 0. (34) 

where n/ {f E u,d, s)^ n,, the particle densities of the quarks and the electrons, respectively. 
The EOS can now be parameterized by only one chemical potential, say fiu = M x. Variation 

omp 

with respect to the trial system Lq gives 

'mF,u,d = mQ- gla = MSu,d (35) 
mF,s =mQ- gl^a = M6s (36) 

In the local equilibrium inside the star the free energy reaches the minimum at a. 

The same result may be achieved calculating the averages of the equation of motions 
(P^) for the effective system Cq. In the mean field approximation the meson field operators 
are replaced by their expectation values. We also consider the isotropic system at rest. As 
< qq >o (calculated with respect to Cq system) depends on the effective quark mass mp (or 
0", a* ) the equation (|Top is highly nonlinear also with respect to a. 

Calculation base on the relation 

dfF ^_ ^ 

=< qq >o 



drriFj 
8 




gives 



FIG. 1: The effective quark masses. 
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The quantum average < g^g/ >o depends on the quark chemical potentials ( pT| , |32D . In the 
result the effective quark effective mass nip also will be dependent on the chemical potentials. 
In the result the solution a of the equation (p^) also will be dependent. The same situation 
will consider other fields. 

The effective quarks mass mpj (or 5f = mpj/M ) dependence on the dimensionless 
Fermi momentum is presented on the Fig.|T]. To calculate the properties of the quark star 
we need the energy-momentum tensor. We define the density of energy and pressure by the 
energy - momentum tensor 

< T^^ >= (P + e)u^,u^ - Pg^y 



where unite vector {u^u^ = 1). Similar to paper UToj [T9| we have introduced the 

dimensionless "Fermi" momentum even at finite temperature which exactly corresponds to 
the Fermi momentum at zero temperature. Both ep and Pp depend on the quark chemical 
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FIG. 2: The quark and neutron (TMl) and hyperon equation of state. 

potential /i or Fermi momentum xp. This parametric dependence on /i (or xp) defines the 
equation of state. The various equations of state for different parameters sets is presented 
on Fig.||. 



The quark star properties 

The numerical results describing the structure of neutron or quark star are based on the 
relativistic mean field theory are now presented. It is possible to describe a static spherical 
star solving the OTV equation. 

dP{r) _ _G P{r) {m{r) + ^^Pjry) 

= Anr'pir) (39) 

dr 

Having solved the OTV equation the pressure p(r), mass m(r) and density p(r) were ob- 
tained. To obtain the total radius R of the star the fulfillment of the condition p{R) = 
is necessary. This allows to determine the total gravitational mass of the star M{R). The 
M{R) for the quark star is presented on the Fig. ^. The R{p) dependence for the quark star 
is presented on the Fig. ||. For the quark star with the central density Pc = 2 10^^ g/crn? the 
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FIG. 3: The M{R) dependence for the quark star. 
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FIG. 4: The R{p) dependence for the quark star. 



star profile in the mean field approach is presented on the Figs. The effective quark 

chemical potential are presented on the Fig. ^ 

The quark and electron Fermi momentum inside the star is presented on Fig. ^. This 
Fermi momentum distribution comes from the quark dispersion relation presented on the 
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1000 




FIG. 5: The effective quark chemical potential for the quark star with the central density pc = 
210^5 g/cm? 

Fig. ^ The meson field and (p^ profile inside the star is presented on the Fig. |^. The 
biggest contributions to the energy density come from the quarks, the gauge boson field u;^, 
0^ and scalar boson field V9, <y9^=. The quark effective mass profile inside the star is presented 
on Fig. ||. This difference between strange and up, down quarks deepens inside the star. Its 
infiuences on the quark partial fraction defined as 

Xf = 1 ^> (40) 

[n^ + + ns) 

where / = {u,d,s), presented on the Fig. We see that both pure neutron star {X^ = 
2/3, Xu = 1/3) and strange star {X^ = X^ = X^) are rather mathematical limits. The 
electron distribution inside the quark star in presented on the Fig. 

In the fiat spacetime the preferred configuration is realized for the system with minimum 
of free energy. As pressure is defined as 

this means that the stable configuration should have bigger pressure. According to the 
equation of state (Fig. H) as was pointed in the strange star configuration is unstable. As 
a consequence there seems to be no chance to find SQM with energies per baryon number 
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FIG. 6: The quark and electron Fermi momentum inside the star 
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FIG. 7: The quark Fermi momentum distribution. 

lower than in ordinary nuclei. This almost rules out the original idea of absolutely stable 
SQM. 

However, when hyperon are included there is crossing point at 



Pe. = 8.92 10^^ ^/cm^ 



13 



100 n 



80 



— 60 
> 
03 

40 

20 









^^^^^ 

Schertler — 


— ^ " • ^ 


Buballa 






0) =8 .w(r) 


' 1 ' 1 


1 


0) =8 „w.(r) 

— n* n,o \ ' 

1 1 1 1 1 



2 4 6 8 10 



r ( km ) 



FIG. 8: The meson field o;^ and (p^ profile inside the quark star. 
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FIG. 9: The quark effective mass profile inside the star. 

above which the strange quark core (the Buballa, Oertel parameters set |Q) may appear. 
This will happen when the star central density pc exceed pcs- However, Fig.H shows that the 
hyperon star is unstable in this region. A mixed hyperon star with strange quark core can 
be stabilized if the the equation of state will be more softer (including the kaon condensation 
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FIG. 10: The quark partial fraction inside the star. 
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FIG. 11: The electron distribution inside the quark star. 



20|| , for example). 

Unfortunately, the strange star with Schetler at al. parameters set is unstable due to 
a large constituent strange quark mass (m* = 549.5 MeV). In the flat Minkowski spacetime 
the nucleon binding energy is presented on Fig.|l2|. 
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FIG. 12: The nuclear binding energy as the relative baryon number u = n/riQ (tiq = 0.17 fm ^ for 
the nuclear symmetric matter). 

All these considerations concerning the star stability ware done in the flat Minkowski 
spacetime. However, a neutron or quark star is mainly binded by gravity. They binding 



energy are equal to 

AE = c^M - c^Mb = {c^M - c^Mp) + {c^Mp - c^Mb), (41) 

where 



- 2Gm{r 
lo ' c^r 

is the proper mass of the compact object and 



Mp = I 4nr'dr{l - P^ y'^p{r) (42) 



M= / A7ir^drp{r) (43) 
Jo 

is the gravitational mass. The gravitational binding energy for the hyperon and quark star 
is presented on Fig. Mb is the baryonic mass deflned as 

Mb = Y.^bNb (44) 

B 

where is the baryon mass and the baryon number is equal to 

Nb = r^Tir^drn - ^^^)-5nB(r) (45) 
JO c^r 

The star mass deficit 

AM = {Mp - M) 
16 




FIG. 13: The star gravitational binding energy as a function of the star baryon number N/Nq 
{Nq = 1.18710^^ is the baryon number of the Sun). 
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FIG. 14: The star mass deficit calculated in the Sun units {Mqc^ = 1.78710^^ ergs) as a function 
of the star central density. 



is presented on Fig.|T^. The quark star is stable with respect to the gravitational binding 
energy (Fig. [l^). The difference of the gravitational binding energy for the hyperon and 
quark star is order of 



8.9 10^2 ergs 



(46) 



Such a doze of energy may be released during conversion of the neutron star into quark 
star. For example, the Mb = Mq hyperon star with the mass M ~ 0.92 Mq, radius 
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R ~ 13km and central density pc ~ Ag/crn? will transform into the quark star with the 
mass M ~ 0.85 M0, radius R ~ 8.5 /cm and central density pc ~ 8.9 g/cm^ (Figs. |13|,|14|). 
This conversion may explain nature of the gamma-ray bursts p2| . 



Conclusion 

Quark meson coupling model was designed to describe both the bulk and internal struc- 
tures of nuclear phenomenology described in terms of the Relativistic Mean Field Theory 
(RMF) The hyperon-hyperon interaction mediated by the (/:>*and (p^ mesons is expected 
to be important for hyperon rich matter present in the cores of neutron star or core of the 
quark star. The quarks interaction is a bit different for the u, d and s quarks. This inter- 
action breaks the starting SU(3) flaver symmetry. The real quark star still is dominated by 
the u, d quarks with smaller presence of the s quarks. Only in the mathematical limit of 
the core high density we have the strange star (X„ = Xd = Xg) case. The appearance of 
deconfined quark matter in the center of the neutron star demands soft nuclear equation of 
state. The quark star occurs stable when the gravitational binding takes into consideration. 

Appendix: The hyperon star in the Relativistic Mean Field Theory 

To compare the quark star with neutron one we shall use the RMF model of the neutron 
star including the A hyperon |^. This is a simple generalization of the paper ||To|. The 
model consists following mesons: ip, u^, with masses = 510 MeV, = 783 MeV, 
Mp = 770 MeV, and v?*, 4>^l with masses m,* = 975 MeV, = 1020 Mel^. The fermion 
fields are composed of protons, neutrons and A hyperon 




with masses M = 938 MeV, Ma 
The Lagrange function is: 



1115.6 Mel/. 



(47) 



(48) 



18 



-\d,ip,d^ip, - \mW. - 1$^,$'- - \m14>^4>^ (49) 
+2^7^D^7/^ - ^(M - g^N^)^ (50) 
+i^^'y''D^^A - V^a(Ma - g^A'f - ga,A^*)i'A (51) 



where F^^ is the stress tensor of the form 



and the tensor _R^j, has the form 



and the scalar potential is 

f/((p) = ^M^^^ + ^g2^^ + -^g^^\ 



where C3, (72, gi are constant of the TMl model [T^. For the hyperon sector 



The covariant derivative are defined now as 



= ("^A* + ^^apNplfT" + ig^N^^p,)ip, (52) 
-D/.V'A = {d^ + ig^A(pf,)i^A (53) 



The coupling constants are 

2 

gaA — '^9ctN 

_V2 

ga*A — gaN 

for SU(6) symmetry, the coupling constants g^rN , gu)N are the same as for the TMl param- 
eters set [|T^ and 

g4>A — — gcuN 

0. The effective mass of nucleon is 

m*,^ = M 5 = M - g„A(T (54) 

and 

m\ = Ma = M 6a = gaA<J - ga,Acr* (55) 

for the A hyperon. The effective mass is presented on Fig |1^. This model allows us to 
calculate the equation of state presented on Fig.^. 
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FIG. 15: The effective mass 6 = m*p/M for nucleoli and A as function of the nucleon Fermi 
momentum x = kp/M. 



[1] J.Boguta, A.R.Bodmer, Nucl. Phys. A292 (1977) 413; Serot B.D., Walecka J.D. ReceM 

Progress in Quantum Hadrodynamics, Int. J. Mod. Phys. E6, (1997) 515-631. 
[2] N.K. Glendenning, Z. Phys., A327 (1987) 295; Also see in Compact Stars by N. K Glendenning, 

( Springer- Verlag, New York, 1997). 
[3] F. Weber Pulsars as Astrophysical Laboratories for Nuclear and Particle Physics, lOP Pub- 

hshing, Philadelphia, (1999). 
[4] A.W.Thomas, Chiral symmetry and the hag model:A new starting point for nuclear physics, 

Adv. Nucl. Phys. 13 (1984) 1-137; P.Guichon, Phys. Lett. B200 (1988) 235; P. Guichon, 



K. Saito, E. Rodionov, A.W. Thomas, The role of nucleon structure in finite nuclei, http 



//xxx. lanl .gov/abs/iiucl-th/9509034| ; H.Shen, H.Toki, Quark mean field model for nuclear 



matter and finite nuclei, [http : //xxx. lanl . gov/ abs/nucl-th/9911046| . 
[5] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 246. 
[6] M. Buballa, M. Oertel,Phys. Lett. B457 (1999) 261; Strange quark matter with dynamically 



20 



generated quark mass, http : / / xxx . lanl . gov/ abs/iiucl-th/98 10529 



[7] K.Schertler, S.Leupold, J.Schaffner-Bielich, Phys.Rev. C60 (1999) 025801; Neutron stars and 



quark phases in the NLJ model, [http : / / xxx . lanl . gov/abs/astro- ph/ 9901152 . 
[8] P. Rehberg, S.P. Klevansky and J. Hiifner, Phys. Rev. C53 (1996) 410. 

[9] I. Zakout, H.J. Jaqaman, S. Pal, H.Stocker, W. Greiner, Hot Hypernudear Matter in the 



Modified Quark Meson Coupling Model, pittp : //xxx . lanl . gov/ abs/iiucl-th/9904084| ; S. Pal, 
M. Hanauske, I. Zakout, H.Stocker, W. Greiner, Neutron star properties in the quark-meson 
coupling model, http : //xxx . lanl . gov/abs/astro-ph/9905010| . 

[10] H. Toki, U. Meyer, A. Faessler, and R. Brockmann, Phys. Rev. C58, 3749 (1998). 

[11] A.R. Bodmer, Phys. Rev. D4 (1971) 1601 

[12] E. Witten, Phys. Rev. D30 (1984) 272. 

[13] E. Farhi and R.L. Jaffe, Phys. Rev. D30 (1984) 2379. 

[14] J. Madsen and P. Haensel (eds.). Strange Quark Matter in Physics and Astrophysics, Nucl. 

Phys. B (Proc. Suppl.) 24B (1991). 
[15] R.Manka, I.Bednarek, J.Przybyla, Phys.Rev C62, 15802, (2000), The neutron star matter in 



a relativistic mean-field theory, [http : //xxx. lanl . gov/ abs/nucl-th/OOOlOlT . 
[16] Y. Sugahara and H. Toki, Prog. Theo. Phys. 92 (1994) 803; H.Heiselberg, Neutron star: recent 

developments, http : //xxx . lanl . gov/abs/nucl-tIi/99120^ (1999). 
[17] A. L. Fetter, J.D.Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, New York, 

(1971). 

[18] R.Feynman, Statistical Mechanics, W.A.Benjamin, Inc. Reading 1972.; R.Manka, J.Kuczyhski, 
G.V.Vittiello, Nucl. Phys. B276 (1986) 533; R.Manka, Annals of Physics, G.V.Vittiello, 199 
(1990) 61. 

[19] Y. Sugahara and H. Toki, Prog. Theo. Phys. 92, 803 (1994). 



[20] N.K. Glendenning, J.Schaffner-Bielich, Phys. Rev. C60 (1999) 025803, [http://xxx.lanl.gov/ 

[ abs/astro-ph/9810290| . 

[21] S. Weinberg, Gravitation and Cosmology, John Wiley & Sons, New York, (1972). 

[22] I. Bombaci, B.Datta, ApJ, 530 (2000) L69. 

[23] N.K. Glendenning, Phys. Lett. B114 (1982) 392; 



21 



